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Resonant control of stochastic spatio-temporal dynamics in a tunnel diode by multiple 

time delayed feedback 

Niels Majer and Eckehard ScholQ 
Institut fiir Theoretische Physik, Technische Universitat Berlin, D-10623 Berlin, Germany 

We study the control of noise-induced spatio-temporal current density patterns in a semiconductor 
nanostructure (double barrier resonant tunnelling diode) by multiple time-delayed feedback. We 
find much more pronounced resonant features of noise-induced oscillations compared to single time 
feedback, rendering the system more sensitive to variations in the delay time r. The coherence of 
noise-induced oscillations measured by the correlation time exhibits sharp resonances as a function 
of r, and can be strongly increased by optimal choices of r. Similarly, the peaks in the power 
spectral density are sharpened. We provide analytical insight into the control mechanism by relating 
the correlation times and mean frequencies of noise-induced breathing oscillations to the stability 
properties of the deterministic stationary current density filaments under the influence of the control 
loop. Moreover, we demonstrate that the use of multiple time delays enlarges the regime in which 
the deterministic dynamical properties of the system are not changed by delay-induced bifurcations. 

PACS numbers: 05.40.-a, 05.45.-a, 72.20.Ht, 72.70.+m 
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I. INTRODUCTION 



It is well known that random fluctuations can seriously 
affect charge transport in semiconductors 1] , which usu- 
ally leads to deterioration of their performance. How- 
ever, recently the constructive role of noise in semicon- 
ductor devices has been recognized. In particular, noise 
was shown to induce coherent radiation in semiconduc- 
tor lasers 0, H, [3] or to induce moving field domains 
in semiconductor superlattices Q whose regularity be- 
comes optimum at some non-zero value of the noise inten- 
sity. This phenomenon is known as coherence resonance 
[li I'^ double-barrier resonant tunneling (DBRT) 
diodes noise can generate spatially inhomogeneous cur- 
rent density patterns in form of breathing current fila- 
ments [8,], however, their regularity decreases monotoni- 
cally with increasing noise intensity, and thus shows no 
coherence resonance. 

The control of the features of noise-induced dynam- 
ics is generally of great importance, and has recently at- 
tracted a lot of attention in the field of nonlinear dy- 
namic systems 0. In [l^ |Tl| a method for manipu- 
lation of essential features of noise-induced oscillations, 
like coherence and time scales, was proposed using a de- 
layed feedback scheme that was originallyused to control 
chaos in purely deterministic systems [l^]- This tech- 
nique was demonstrated to be effective for control of 
noise-induced oscillations in either simple (generic) sys- 
tems [isl nil. Iisl HI. IitI ITsI or more complex, spatially ex- 
tended systems [1911201 [2l|. For the DBRT nanostructure 
it was shown that time-delayed feedback can either in- 
crease of decrease the regularity of noise-induced breath- 
ing filaments and, moreover, can even lead to spatial ho- 
mogenization of current density patterns (2^ . For de- 
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terministic systems the original single time delayed feed- 
back was extended by using multiples of the delay time r 
weighted with a memory parameter [23| , which generally 
leads to lar ger control domains and more efficient con- 
trol [H, IHillE I13- ^ simple stochastic Van der Pol 
oscillator this has also been shown to yield drastically 
increased correlation times psl j. 

In the present work we study the effect of multiple time 
delayed feedback on the stochastic spatio-temporal pat- 
tern formation in the DBRT model. Compared to the 
single time delayed feedback in the same system [22j , we 
find much sharper resonances of the spectral and cor- 
relation properties in dependence upon the delay time. 
For parameter values close to, but below, a Hopf bifurca- 
tion we show that in these sharp, pronounced resonances 
the temporal regularity is significantly increased and the 
power spectral width becomes much narrower. Moreover, 
we demonstrate that the use of multiple time delays en- 
larges the control parameter regime in which the original 
deterministic dynamical properties of the system do not 
change, i.e. the delay-induced bifurcations occur only at 
larger feedback strength. 

The paper is organized as follows. In Section II the 
DRBT model is described, and the dynamical properties 
of the system for our chosen parameters are discussed. 
Section III is devoted to the effects of the multiple time 
delayed feedback upon noise-induced dynamics, and in 
Section IV we draw conclusions. 



II. MODEL 

In our study, we use a deterministic model for the 
DBRT suggested in [29| and add two sources of random 
fiuctuations as proposed in Q. Furthermore we use the 
time delayed feedback scheme which was already applied 
to this system in [23 | and extend it in order to take mul- 
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tiple time intervals into account in the feedback loop: 



da{x, t) 
dt 

du{t) 



f{a,u) 



(1) 



^-{Uo-u- rJ) + DuTl{t) + F{t) 



where all quantities are dimensionless. The dynamical 
variable a{x, t) describes the charge carrier density in- 
side the quantum well, whereas u{t) is the voltage drop 
across the device. The spatial coordinate x denotes 
the direction perpendicular to the current flow, and t 
is time. In terms of nonlinear dynamics, in the determin- 
istic {Da — Du = 0) and uncontrolled {K = 0) case, this 
is a reaction-diffusion model of activator-inhibitor type, 
where a is the activator, and u is the inhibitor [soj . 

The net tunneling rate of the electrons through the 
two energy barriers into and out of the quantum well is 
modelled by the nonlinear function [29j : 



f{a,u) 

Jin ' 
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(2) 



where d is the effective thickness of the double-barrier 
structure, re = {Aneeofi^) /{e'^m) is the effective Bohr 
radius in the semiconductor material, e and eg are the 
relative and absolute permittivity of the material, and 
Xq and 7 describe the energy level and the broadening 
of the electron states in the quantum well and rj,, is the 
dimensionless Fermi level in the emitter, all in units of 
kBT. Throughout the paper we use values of 7 = 6, 
d/vB — 2, rje — 28 and xq = 114, corresponding to typical 
device parameters at 4 K ^29*]. 

The effective diffusion coefficient D{a) resulting from 
the inhomogeneous lateral redistribution of carriers and 
from the change in the local potential due to the charge 
accumulated in the quantum well by Poisson's equation 
is given by [3]| : 



Dia)^a[^ + - ^-j— 

\rB l-exp(-aj 



(3) 



It describes the diffusion of the electrons within the 
quantum well perpendicular to the current flow. J = 
T lo ^'^^ gives the total current through the device, 
where j(a,u) = ^ {jt„ + jout) = 5 (/(a, u) -I- 2a) is the 
local current density within the well. The system's width 
is fixed at a value of L = 30 and homogeneous Neumann 
boundary conditions are used. 

In Eq. IT]) we use uncorrelated Gaussian white noise 
sources ^(2;, t) and ri{t) with noise intensities Da and Du- 



(C(x,t)) = (,y(t))=0 

{a^,t)ax',t'))^6ix 

{v{t)v{t')) = S{t - 



ixe[0,L]), 
~x')5{t-t'), 
t'). 



(4) 



Physically, D^ can be realized by an external tunable 
noise voltage source in parallel with the supply bias. Da 
describes internal fluctuations of the local current density 
which could be caused, e. g., by shot noise [l|. 

The control force F(t) represents a control voltage 
which is constructed recursively from a time-delayed 
feedback loop with delay time r, feedback strength K > 
0, and memory parameter R, and can be written as 

F{t) = K{u{t-T) -u{t))+RF(t~T) (5) 
^ K^R''[u{t-{n+l)T)~u{t-nT)]. (6) 

n=0 

The first Eq. ^ is the local balance equation of the 
charge in the quantum well, and the second equation rep- 
resents Kirchhoff 's law of the circuit in which the device 
is operated. The control parameters are the external bias 
voltage Uq, the dimensionless load resistance r, and the 
time-scale ratio e = RC/Ta, which is related to the load 
resistance R, and the parallel capacitance C of the at- 
tached circuit, normalized by the tunneling time Tq. A 
discussion of the various deterministic bifurcation scenar- 
ios can be found in [2^, [2^ [sl] . 

We fix e = 6.2 slightly below the Hopf bifurcation, 
which occurs at £Hopf ~ 6.469. In this regime we have 
two fixed points: (i) a stable, spatially inhomogeneous 
fixed point and (ii) a spatially homogeneous fixed point 
which is stable with respect to completely homogeneous 
perturbations but generally unstable against spatially in- 
homogenous fluctuations (saddle-point). Although the 
deterministic systems rests in the inhomogeneous steady 
state, noise can induce irregular spatio-temporal oscilla- 
tions of the current density [1] . In the following we shall 
study how these noise-induced oscillations are influenced 
by the control force. 



III. MULTIPLE TIME DELAYED FEEDBACK 
CONTROL 

Fig[T]shows simulations of the spatio-temporal dynam- 
ics under the influence of noise and delayed feedback. 
The voltage time series (a), the spatio-temporal charge 
density patterns (b), and the current-voltage projection 
of the infinite-dimensional phase space (c) are depicted. 
Noise induces small spatially inhomogeneous oscillations 
around the inhomogeneous steady state (breathing cur- 
rent filaments). In the J-u phase portrait (c), the spa- 
tially inhomogeneous steady state (fixed point) is deter- 
mined by the intersection of the load line (null isocline 
ii — Q, blue dash-dotted) with the nuUcline d = for 
inhomogeneous a{x,t) (red dotted). The neighboring in- 
tersection of the load line with the nuUcline d = for ho- 
mogeneous a (black solid curve) defines the second, spa- 
tially homogeneous fixed point which is a saddle. With 
increasing noise intensity (Fig|2|) the oscillation ampli- 
tude becomes larger, the oscillations become more irreg- 
ular, and finally, at even larger noise, the oscillations are 
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more spatially homogeneous, i.e., in the phase space they 
are more centered around the homogeneous fixed point 
(FigEl). 
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FIG. 1: (color online) Stochastic spatio-temporal dynamics 
under multiple time-delayed feedback control, (a) Voltage 
time series u(t) (in units of 0.35 mV), (b) charge carrier den- 
sity a{x, t) (in units of 10^''/cm'^), (c) phase portrait of current 
J (in units of 500 A/cm^) vs. voltage u. Space x and time t 
are scaled in units of 100 nm and 3.3 ps, respectively, corre- 
sponding to typical device parameters at 4 K 29] . Parameters 
are Uo = -84.2895, r = -35, e = 6.2, D„ = 0.1, Da = 10"'', 
K = Q.l, T = 6.3, ii = 0.5. 
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FIG. 2: (color online) Same as Fig.[T]for Du = 1.0. 

We shall now investigate how the regularity and the 
time-scales of these noise-induced oscillations depend 
upon the feedback parameters if, i?, r. 



A. Linear stability analysis 

To this purpose we first examine the stability prop- 
erties of the inhomogeneous fixed point (ao(a;),uo) un- 
der the influence of the control force. We perform a 
linearization of the original continuous system ([1]) (with 
Du = Da = 0) around the inhomogeneous fixed point 
along the same lines as in [l^, [3^ . We use an exponential 
ansatz for the deviations from the fixed point 5a{x^ t) 



a(x, t) — ao(x) 



d{x) and Su{t) = u{t) — uq 
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FIG. 3: (color online) Same as Fig.[T]for Du = 2.0. 

The resulting coupled eigenvalue problem reads 
Aa{x) — Ha{x) + fu{x)u, 



Au = 



eL 



K- 



ja{x)a{x)dx i — -u 



-At 



1 -Rc 



-A: 



(7) 



(8) 



where we have introduced a self-adjoint linear operator 
H. Its eigenfunctions ^'^ and eigenvalues correspond 
to the voltage-clamped system, Su = 0. Furthermore, 



fu 



du 

Ju — 



]a 



dj_ 

da 



(9) 



1 



•dj_ 
L Jq du 



dx. 



Due to the global constraint, Eq.® mixes the eigen- 
modes 4'^ and both equations have to be solved simul- 
taneously. An expansion of the eigenmodes a of the full 
system in terms of the eigenmodes of the voltage- 
clamped system, keeping only the dominant eigenmode 
with eigenvalue Aq > 0, leads to a characteristic equa- 
tion for the eigenvalues A 



A^ 



l + rJu 



-Ao A 



+ (A-Ao)if^ 



-At 



Ao, 



1 - Rc- 



__^(l + ^^,) = 0, (10) 



where the static differential conductance at the inhomo- 
geneous fixed point ad = ^\ has been introduced. 
In [24I a more detailed derivation of the characteristic 
equation is given for the special case R~Q. The exten- 
sion to the case i? 7^ is straightforward. 

Without control, K = eq. pO|) reduces to a char- 
acteristic polynomial of second order, which gives the 
well-known conditions for stability of a filament (ssj 



A^i±^-Ao>0, 



C 



e 
Ao 



(11) 



{l + rad) > 0. 
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and a Hopf bifurcation occurs on the two-dimensional 
center manifold if ^ = 0. 

With control, eq. (jlOp can be expressed as 

+ Ak + {K- B)K-^—^—^ + C (12) 

with S = Aq > 0. The parameters A, B, C can be 
calculated directly from dfl]) yielding A = 0.0447, 
B = 1.0281 and C = 1.1458. 

Using Eq. (fT^ we can calculate the domains of stabil- 
ity in the t-K plane numerically for selected values of the 
memory parameter R. In order to find the curves con- 
taining the boundaries of stability of the inhomogeneous 



fixed point as a subset, we set A = p + iq with p = and 
separate Eq. (jl2p into real and imaginary parts: 

[BK -R{C - q^)] cos(gT) - q [AR + K) sin(gr) 

= BK + {q^ - C) (13) 

q {AR + K) cos(gT) + [BK + R {q^ - C)] sin(gT) 

= qiA + K) (14) 

Using Eq. dTS]) and Eq. ([14]) the boundary of stabil- 
ity can be obtained from the set of parametric functions 
K{q) and T{q) using q = /m(A) as the curve parameter. 



^ (AV + (c-g^)^)(i + i?) 

2BC -2{A + B)q^ 

, , 1 / . / q{AB + C - q'){l - R)K \ \ 

T{q) - - (^arcsm ^ _ ^2)2)^2 + 2{-BC +{A + B)q^)RK + [B^ + q^)K^ ) + j 
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FIG. 4: (color online) (a)-(d) Stability domains of the inho- 
mogeneous fixed point in the t-K plane of the deterministic 
system ([1} {Du = Da =0), for selected values of the memory 
parameter R. Blue lines: Solutions of p2|l with f?e(A) = 
calculcated from Eq. (|15|l . Orange (shaded) region: regime 
of stability of the fixed point obtained from the numerical 
solution of Eq. ([l]). Black horizontal line: upper bound for 
K where the fixed point is stable for all values of r, calcu- 
lated from Ea. p6|l . The black diamond in panel (c) marks 
the parameter values for which a J — u phase portrait of the 
delay-induced limit cycle is shown in the inset. 



branches. The boundaries are in good agreement with 
the domain of stability obtained from dynamical simula- 
tions of the nonlinear system equations IH) {Du ~ Da = 
0), shown as orange (dark shaded) areas. The inset of 
panel (c) shows the delay- induced limit cycle in the J — u 
phase space for parameters outside the stability domain 
of the inhomogeneous fixed point. The stability domains 
increase significantly with increasing memory parameter 
R from (a) to (d) . The modulation of their boundaries in 
dependence on r results from the cross-over of differerent 
eigenvalue branches, which is a typical feature of delay 
differential equations. 

From Eq. (|15p it is possible to calculate an upper 
bound Kc of K for which the stability properties of 
the uncontrolled deterministic system remain unchanged 
over the whole r interval, meaning that no delay-induced 
Hopf bifurcation occurs. 



Kr = 



A^{A + B) [BC -G) + {AC + Cf 

2{A + BfG 
0.1059(1 + i?) 



(1 + i?) 



G: = ^/A^G{B{A + B) + C) 



(16) 



Fig. [4] shows this upper bound plotted as black horizontal 
line. 



Fig. 0] shows these curves, Eq. p^ . as blue lines. The 
boundaries of stability, where a delay-induced Hopf bi- 
furcation of deterministic breathing oscillations occurs, 
are a subset of these curves, because the fixed point may 
already be unstable when a complex eigenvalue crosses 
the imaginary axis, due to other unstable eigenvalue 



B. Correlation times 

To quantify the temporal regularity of the noise- 
induced oscillations, we evaluate the correlation time [341 
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FIG. 5: (color online) (a)-(d) Power spectral density Suu{f) 
of the dynamical variable u in dependence of the frequency 
/ and the delay time t for selected values of R {K = 0.1, 
e = 6.2, Du = 0.1, Da = 10"*). 

calculated from the voltage signal, 

tcor = \ I (17) 
Jo 

I 



which is shown as black curves in Fig. [6l 

At certain resonant values of r (left column) the spec- 
tral peaks become extremely sharp. With increasing 
memory parameter R the broad spectra prevail over 
larger intervals of r, whereas the regime of sharply 
peaked spectra becomes smaller. Thus multiple time 
feedback control exhibits much more pronounced reso- 
nant features both in the frequency domain and in the 
delay time. Since a sharply peaked spectrum gives rise 
to long correlation times (which are, in the linear regime, 
proportional to the inverse spectral width) we expect 
the domains of strong correlation to shrink with increas- 
ing memory parameter and the domains of low corre- 
lation to increase. Extracting from the Fourier power 
spectral density the autocorrelation function ^'(s) = 
!'^^Suu{f)e^''^f' df and using Eq. ^ we obtain the 
correlation time tcor in dependence of r. This is shown 
in Fig. [7] for different values of the memory parameter 
R. The feedback strength is kept at a constant value of 
K = 0.1, where the system is below the Hopf bifurcation 
for all values of t and R. 

For small memory parameter R the correlation times 



where \l/(s) = {{u{t) — {u)) {u{t + s) — {u)))^ is the au- 
tocorrelation function of the variable u{t) and — ^'(0) 
its variance. 



In order to investigate the influence of multiple time 
delayed feedback control in the DBRT, we systemati- 
cally study the dependence of the correlation time from 
Eq. p7|) upon the control force parameters r, if, and R. 



In Fig. [5] the Fourier power spectral density Suu ob- 
tained from the time series u{t) is shown in dependence 
of the delay time t for different values of the memory pa- 
rameter R. The shape of the spectra S'„„(27r/) alternates 
between broad and sharply peaked with varying r. This 
shows up more clearly in the corresponding sections at 
fixed T depicted in Fig. [S) A very good analytic approx- 
imation of the power spectral density can be obtained 
by a straightforward extension of the argument in [25| to 
multiple time-delayed feedback: 



(18) 



alternate between high and low values with growing r. 
For higher memory parameters R the peaks in correla- 
tion time indeed become narrower and sharpen up, and 
the domains of low correlation time increase. The sta- 
bility of the inhomogeneous fixed point reveals a rela- 
tion between properties of the controlled deterministic 
system and the noise-induced dynamics: maximum reg- 
ularity of noise-induced oscillations is attained when the 
deterministic fixed point is least stable. This feature is 
maintained for all values of the memory parameter R. 
In the case of small R the crossover of the real part of 
eigenvalue branches also determines the location of the 
minima in correlation time. In that case two eigenmodes 
with the same stability (real part) but different frequen- 
cies are present in the system, resulting in rather irreg- 
ular noise- induced dynamics. For large memory param- 
eters the broad domains of low correlation display many 
eigenmodes that are not well separated (stability-wise) 
causing irregular mixed dynamics. 

The regularity of the noise-induced oscillations in de- 
pendence of the control force strength K is visualized in 
Fig. [5] The correlation time vs K is shown for different 
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FIG. 6: (color online) Power spectral density Suuif) of the 
dynamical variable u in dependence of the frequency / for 
various delay times r and memory parameters R {K — 0.1, 
€ = 6.2, = 0.1, Da = 10"*). 



values of R and r. Depending upon the chosen value of r 
the correlation time increases or decreases with growing 
K. An optimal value of r leads to more regular oscilla- 
tions whereas a non-optimal value of r results in more 
irregular oscillations. In the case of an optimally cho- 
sen r and a value oi K > 0.1 the curves split for larger 
K, and the one with largest R (blue full squares) attains 
the highest correlation times, clearly above the curves for 
R = (red plus) and R = 0.5 (violet asterisks). Compar- 
ing the values of tcor at the threshold of delay-induced 
Hopf bifurcation of the fixed point, marked by K\j^pp 
the order of increasing correlation time is from small to 
large R. 

The dependence of the correlation time on the mem- 
ory parameter R is shown in Fig. [H^a) at a fixed value 
of T = 6.19 and values of K — 0.1 and 0.2. Initially 
starting at a higher value of tcor the correlation time in- 
creases more rapidly and reaches higher values in the case 
of larger K (K=0.2). Figure [51[b) shows the domain of 
stability in the K-R plane. The colored (dark) region 
is the domain where the inhomogeneous fixed point is 
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FIG. 7: (color online) Correlation times tcor (upper panels) 
and deterministic stability of the inhomogeneous fixed point, 
i?e(A) (lower panels), in dependence of the delay time r for 
different values of the memory parameter R (a)-(d) and fixed 
K = 0.1. The red (dark) curves in the lower panels mark the 
leading eigenvalue, which governs the overall stability of the 
fixed point. Parameters: e = 6.2, Du = 0.1 and Da = 10~* 
(in the panels showing tcor)- 
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FIG. 8: (color online) Correlation times tcor in dependence of 
K for different values of the memory parameter R and optimal 
and non-optimal r {e — 6.2, Du = 0.1, Da — 10"''). Kjfc,pf, 
Kjjapf a-nd Kfjc,pf mark the values of K at which the Hopf 
bifurcation occurs for R=0, R=0.5, and R=0.9, respectively 
(e = 6.2, Da = 0.1, Da = 10"*). 



stable under the influence of the feedback loop without 
noise. The two lines mark the chosen values of K used 
in Fig. [DJa). Again, the correlation time is maximum if 
the fixed point is closest to its boundary of stability. 



C. Mean period 

Here we will establish a relation between the linear 
modes of the inhomogeneous fixed point and the time 
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FIG. 9: (color online) (a): Correlation times tcor in depen- 
dence of the memory parameter R vor different values of the 
coupling strength K and r = 6.19. (b): Domain of stability of 
the inhomogeneous fixed point (green shaded) in the {K, R) 
plane for r = 6.19. The horizontal lines mark the values of K 
chosen in (a) (e = 6.2, L»„ = 0.1, Da = lO"*). 

scales of the noise-induced oscillations in its vicinity. 
The stability and angular frequency of the eigenmodes 
are given by the real and imaginary parts of the eigen- 
values Ki of Eq. (fT2|l . Re{Ki) and Im(Ai), respectively. 
Fig. (fTO|l shows a plot of the main periods of the noise- 
induced oscillations obtained from the main peak in the 
Fourier power spectral density (blue open circles) com- 
pared to the periods calculated from the imaginary part 
Im(Ai) of the eigenvalues of the deterministic fixed point 
using Eq. (fT2|) . The red (solid) curve shows the pe- 
riods Tq = 27r//m(Ao) of the leading eigenvalue. For 
small memory parameter R the period of the oscilla- 
tions closely follows the eigenperiod corresponding to the 
leading eigenvalue over the whole r interval, whereas for 
larger R this feature is only maintained in a narrow do- 
main where the leading eigenvalue has a real part close to 
zero and is clearly separated from all other eigenmodes 
(compare Fig. [T]). 
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FIG. 10: (color online) Main period To = 1/fmax of the noise 
induced oscillations obtained from the main peak in the power 
spectral density (blue open circles). Red (solid) curve: Pe- 
riod T = 27r//m(Ao) of the leading eigenvalue. Green (grey) 
curves: Period T = 27r/7m(Ai) (K = 0.1, e = 6.2, Du = 0.1, 
Da = 10-*). 

IV. CONCLUSION 

In conclusion, we have shown that multiple time- 
delayed feedback control leads to more pronounced res- 
onant features of noise-induced spatio-temporal current 
oscillations in a semiconductor nanostructure compared 
to single time feedback. The regularity of noise-induced 
oscillations measured by the correlation time exhibits 
sharp resonances as a function of the delay time t, and 
can be strongly increased by control with optimal choices 
of r, whereas it decreases in a broad range of non-optimal 
values of t. Thus the system is more sensitive to vari- 
ations in T. Similarly, the peaks in the power spectral 
density are sharper and exhibit stronger resonances in de- 
pendence on T for multiple time feedback, whereas their 
position, i.e., the main period of the oscillations, is less 
sensitive to variations in r in wider intervals. 

The regularity and time scales of noise-induced breath- 
ing oscillations are related to the stability properties of 
the deterministic stationary filamentary current pattern 
(fixed point) under the infiuence of the control loop. 
Maximum regularity is attained if the fixed point is least 
stable in the deterministic case. In the domains of high 
temporal correlation the period of the noise-induced os- 
cillations corresponds to the eigenperiod of the leading 
eigenvalue in the deterministic system. 

Furthermore, we have shown that using multiple time 
delayed feedback control compared to single time-delayed 
feedback leads to larger regimes of stability of the station- 
ary filamentary current pattern in the deterministic sys- 
tem, and delay-induced bifurcations occur only at larger 
values of the control amplitude K. 
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